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ABSTRACT 

The Leo Group elliptical NGC 3379 is one of the few normal ellipti- 
cal galaxies close enough to make possible observations of resolved stellar 
populations, deep globular cluster (GC) photometry and high signal-to-noise 
GC spectra. We have obtained Gemini/GMOS spectra for 22 GCs associated 
with NGC 3379. We derive ages, metallicities and a-element abundance ratios 
from simple stellar population models using the multi-index \ 2 minimisation 
method of Proctor & Sansom (2002). All of these GCs are found to be con- 
sistent with old ages, i.e. > 10 Gyr, with a wide range of metallicities. This is 
comparable to the ages and metallicities Gregg et al. (2004) find for resolved 
stellar populations in the outer regions of this elliptical. A trend of decreasing 
a-element abundance ratio with increasing metallicity is indicated. 

The projected velocity dispersion of the GC system is consistent with 
being constant with radius. Non-parametric, isotropic models require a sig- 
nificant increase in the mass-to-light ratio at large radii. This result is in 
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contrast to that of Romanowsky et al. (2003) who find a decrease in the ve- 
locity dispersion profile as determined from planetary nebulae. Our constant 
dispersion requires a normal sized dark halo, although without anisotropic 
models we cannot rigorously determine the dark halo mass. 

A two-sided \ 2 test over all radii, gives a 2a difference between the mass 
profile derived from our GCs compared to the PN-derived mass model of 
Romanowsky et al. (2003). However, if we restrict our analysis to radii be- 
yond one effective radius and test if the GC velocity dispersion is consistently 
higher, we determine a > 3cr difference between the mass models, and hence 
favor the conclusion that NGC 3379 does indeed have dark matter at large 
radii in its halo. 

Key words: globular clusters: general - galaxies: individual: NGC 3379 - 
galaxies: star clusters - galaxies: kinematics and dynamics. 

1 INTRODUCTION 

The Globular Cluster (GC) systems of numerous elliptical galaxies have been well-studied 
photometrically (e.g. Gebhardt & Kissler-Patig 1999; Kundu & Whitmore 2001; Larsen et 
al. 2001). Since the discovery of the bimodal colour distribution of GCs in elliptical galaxies 
(see Harris 1999 for a review), several models have been proposed to explain how these GC 
systems are formed. 

Ashman & Zepf (1992) propose that some, or all, of the high- met allicity GCs will be 
formed during gas-rich merger events and therefore be of a similar age to those mergers. For 
the case of a series of minor gaseous mergers, a roughly monotonic increase in metallicity is 
expected with decreasing age. 

The multi-phase collapse model of Forbes, Brodie & Grillmair (1997) proposes that the 
majority of GCs are native to their galaxy and formed during two, or more, proto-galactic 
collapse phases. Red, metal-rich, GCs are expected to be 2-4 Gyrs younger than their blue, 
metal-poor counterparts according to the semi-analytic GC formation model of Beasley et 
al. (2002). 

Cote, Marzke & West (1998), invoke tidal capture of GCs from dwarf galaxies to account 
for the blue GC population of large elliptical galaxies. A higher fraction of blue, metal-poor, 
GCs is expected compared with spirals of a similar luminosity. Assuming that blue GCs 
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form at the same time as their parent galaxy, "down-sizing" (Cowie et al. 1996; Kodama et 
al. 2005) suggests that blue GCs captured from dwarf galaxies should be younger than the 
native blue GCs. 

Thus a key discriminant of the formation models is GC age. Unfortunately from pho- 
tometry alone it is very difficult to determine the age and metallicity of individual GCs, due 
to the age-metallicity degeneracy (Brodie et al. 2005), therefore high quality spectra are re- 
quired. The stellar population properties of individual GCs can be measured accurately with 
low-resolution integrated spectra of sufficient signal-to-noise (>30). The age-metallicity de- 
generacy can be broken by measuring spectral Lick indices (Worthey 1994), and comparing 
those indices to single stellar population (SSP) models. High signal-to-noise spectroscopic 
studies (i.e. those capable of determining ages for individual GCs) of elliptical galaxy GC 
systems are limited to around a dozen galaxies (e.g. NGC 1399 Forbes et al. 2001; NGC 1023 
Larsen & Brodie 2002; NGC 1316 Goudfrooij et al. 2001; NGC 3610 Strader et al. 2003, 2004; 
NGC 2434, NGC 3379, NGC 3585, NGC 5846 and NGC 7192 Puzia et al. 2004; NGC 5128 
Peng, Ford & Freeman 2004; NGC 4365 Brodie et al. 2005; NGC 1052 Pierce et al. 2005). 
Except for Goudfrooij et al. and Peng et al. , these works find the majority of GCs to be 
old (>10 Gyrs) with a small fraction of young GCs for some galaxies. 

In this work we report the analysis of Gemini/GMOS spectra of a sample of GCs in 
NGC 3379. This is a typical, nearby (10.8±0.6 Mpc, Gregg et al. 2004) E0/1 elliptical 
galaxy. It is of intermediate luminosity (My=-21.06) with typical early-type colours, Mg2 
index and velocity dispersion (Davies et al. 1987) for its luminosity. There is no sign of any 
optical disturbance (Schweizer & Seitzer 1992). There is, however, a small nuclear dust ring 
at a radius of 1.5" (van Dokkum & Franx 1995) and some ionized gas that extends to a 
radius of 8" (Macchetto et al. 1996). 

NGC 3379 is one of only a few elliptical galaxies to lie close enough that its stars can be 
resolved by the Hubble Space Telescope (HST). Resolved stellar population measurements 
using the iifSTNICMOS camera in the Jand # bands have been made by Gregg et al. (2004). 
Measuring individual stellar magnitudes and colours to just below the Red Giant Branch 
(RGB) tip, they found the outer region stellar population to be old, with ages > 8 Gyr 
and a mean metallicity around solar. They noted similarities to the Galactic bulge. For the 
central region, Terlevich & Forbes (2002) used Lick indices to estimate an age of 9.3 Gyrs, 
[Fe/H]=+0.16 and [Mg/Fe]=+0.24. This suggests that radial age and metallicity variations 
are relatively small. 



4 Pierce et al. 

The U, B, R, I photometry of the GC system of this galaxy reveals the classic GC 
colour bimodality (Whitlock, Forbes & Beasley 2003). Whitlock et al. estimate a GC specific 
frequency of Sat = 1.1±0.6 which is low for an elliptical galaxy, although it is important to 
remember that NGC 3379 belongs to a relatively small group of galaxies (Leo) and therefore, 
a low S N is not unexpected (Bridges 1992). The B, V,R photometry of Rhode & Zepf (2004) 
indicates a similarly low S N of 1.2±0.3 and that approximately 70% of NGC 3379 GCs are 
blue. 

To derive GC ages and metallicities we measure Lick indices from our GMOS spectra. 
We then apply the multi-index ^-minimisation method of Proctor & Sansom (2002), which 
employs all the available Lick indices to break the age-metallicity degeneracy and simulta- 
neously measure the a-element abundance ratios (see Proctor, Forbes & Beasley 2004 for 
an application to Galactic GCs and Pierce et al. 2005 for NGC 1052 GCs). 

In addition to the measurement of line indices, integrated spectra of GCs allow the 
measurement of recession velocities, which confirm whether the GCs are indeed associated 
with the galaxy being studied. With a large enough sample of measured GC velocities it is 
possible to probe the gravitational potential and therefore the dark matter halo of a galaxy 
(e.g. Zepf et al. 2000; Cote et al. 2001; Cote et al. 2003; Peng et al. 2004; Richtler et al. 
2004). 

Measurements of planetary nebulae (PN) kinematics in NGC 3379 (Romanowsky et 
al. 2003) reveal a decreasing velocity dispersion profile at large radii. This suggested that 
NGC 3379 has a significantly lower mass-to-light ratio (M/L = 7.1±0.6 at 5 R e ) than most 
large ellipticals, indicating a minimal Dark Matter (DM) halo. This result would appear 
to be in conflict with the standard Cold Dark Matter (CDM) galaxy formation picture in 
which all galaxies lie within significant DM halos. Here we use the measured velocities of 
our sample of GCs, as well as the GC data of Puzia et al. (2004), as an alternative data set 
to test this important claim. 

In this paper we present our GMOS observations and data reduction of GC spectra in 
Section 2. An analysis of ages, metallicities and a-element abundance ratios derived from 
Lick indices is presented in Section 3. In Section 4 we focus on GC kinematics and the 
implications for the dark matter halo of NGC 3379. We discuss our results in Section 5 and 
present our conclusions in Section 6. 



Globular Clusters in NGC 3319 5 
2 OBSERVATIONS AND DATA REDUCTION 



The observations described below are part of Gemini program GN-2003A-Q-22. GC can- 
didates were selected from Gemini North Multi-Object Spectrograph (GMOS; Hook et al. 
2002) imaging, obtained during January 2003, for three fields around NGC 3379. The data 
reduction and GC candidate selection process is almost identical to that of Forbes et al. 
(2004) for NGC N4649, briefly, SExtractor (Bertin & Arnouts 1996) was used to select GC 
candidates based upon their Gemini zero point colours (0.7< g — i <l-4), magnitudes (i < 
24) and structural properties (i.e. objects with stellarity index >0.35). Interstellar reddening 
towards NGC 3379 is E(B-V)=0.024 mag and is not taken into consideration (Schlegel et al. 
1998). Spectroscopy of NGC 3379 globular clusters were obtained with GMOS on the Gem- 
ini North telescope in the months of February and April 2003. GMOS masks for three fields 
were designed, but only the central field was observed within the time allocation. Seeing was 
typically ~0.9 arcsec. 

The GMOS CCDs consist of 3 abutted 2048x4608 EEV chips with a plate scale of 
0.0727 arcsec/pixel (un-binned). For our set-up, we binned 4x in the dispersion direction, 
yielding 1.84A per binned pixel, giving a final resolution of FWHM ~4A. The dispersion runs 
across the detectors in the GMOS instrument, resulting in two ~20A gaps in the dispersion 
direction of the spectra. The B600_G5303 grating was used, with central wavelengths for 
successive exposures set alternately to 5000A and 5025A (in order to obtain full wavelength 
information across the gaps in the GMOS detectors). The effective wavelength range of each 
slit-let is a function of its position on the mask, but typically covered 3800A-6660A. A slit 
width of 1.0 arcsec was chosen to match the seeing and the minimum slit length was 6 
arcsec, a compromise between maximising the number of slit-lets and allowing for reliable 
sky-subtraction. Exposures of 20 x 1800s were taken, yielding a total of 10 hours on-source 
integration time. Bias frames, dome flat-fields and Copper- Argon (CuAr) arc exposures were 
taken as part of the Gemini baseline calibrations. 

These data were reduced using the Gemini/GMOS packages in IRAF and a number of 
custom made scripts. From the CuAr arcs, wavelength solutions with typical residuals of 
O.lA were achieved. Objects and sky regions in the object spectra were manually identified 
in cross-sections of the two-dimensional images. The spectra were then extracted by tracing 
them in the dispersion direction. After some experimentation, optimal (variance) extraction 
was found to yield the best results since our data are over-sampled on the detector. In 
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some cases, objects were too faint to trace individually and we therefore co-added several 
2-d images, taken adjacent in time, to act as a reference for the extractions. We verified 
beforehand that flexure was minimal between the reference images. Finally, the extracted 
spectra were median combined and weighted by their fluxes with cosmic ray rejection. 

In the absence of any velocity standard stars, the recession velocities were measured by 
using four Bruzual & Chariot (2003) model stellar energy distributions (SEDs) for old ages 
with metallicities [Fe/H] = -2.25, -1.64, -0.64 and -0.33. The task FXCOR in IRAF was 
used and the average (weighted by the cross-correlation peak-height) was taken. Objects with 
recession velocities in the range 900±400 km/s are potentially associated with NGC 3379. 
These are presented in Table 1 and also used in our kinematic analysis in Section 4. There 
was one foreground star and one background galaxy out of the 24 spectra obtained. Our 
low contamination rate of 9% is due to good imaging and colour selection. We note that 
the spectrum of gl420 shows emission lines, most notable 4959 and 5007 A [OIII] which we 
assume are due to a planetary nebula (see Fig. 1). Minniti & Rejkuba (2002) report a similar 
discovery for a GC in NGC 5128 and provide a discussion of the implications. Another GC, 
gl426, shows unexplained emission features around ~5000 A. 

A first-order flux calibration was carried out by normalising the spectra with a low-order 
polynomial. In order to measure Lick indices, we convolved our spectra with a wavelength- 
dependent Gaussian kernel to match the resolution of the Lick/IDS system (see Beasley et al. 
2004b). Lick indices (Trager et al. 1998) were measured from our normalised spectra. Due to 
the variable wavelength ranges in these spectra, the same set of indices could not be measured 
for all spectra. However, all covered a wavelength range of 4500-5500 A. Uncertainties in 
the indices were derived from the photon noise in the unfluxed spectra. No Lick standard 
stars were observed so we therefore cannot fully calibrate the GCs onto the Lick system. 
Consequently there are some systematic differences between some of the measured indices 
and those used in the SSP models. These issues are discussed further in Sections 3 and 5. 
Measured line indices and uncertainties are presented in Tables 2 and 3. 

The final spectra have S/N = 18-58 A -1 at 5000 A, giving errors in the FL3 index of 
0.13-0.44 A (we note that gl566 has an H/3 error of ±0.065 A and S/N=118, however this 
is not representative as this GC is ~1.5 mags brighter than any other). 
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Figure 1. Normalised GC spectra which have been offset by one unit. These spectra have not been de-redshifted. These sample 
spectra show the wavelength range that the majority of our spectra cover and display the range of S/N and metallicity present. 
The spectrum of gl426 shows emission features around 5000 A and 5890A of unknown origin. The PN hosting GC, gl420, is 
also plotted and the emission due to 4959 and 5007 A [OIII] lines can be seen redshifted to 4971 and 5019 A. 

3 AGES, METALLICITIES AND a-ELEMENT ABUNDANCE RATIOS 

In this section we describe how we derive ages, metallicities and a-element abundances. The 
resulting values are listed in Table 4. 

We apply the \ 2 multi-index fitting technique of Proctor & Sansom (2002) for this 
analysis. This method involves the comparison of the measured Lick indices with SSP models 
(its application to extra-galactic GCs is described fully in Pierce et al. 2005). The SSP models 
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Table 1. Confirmed Globular Clusters around NGC 3379. Cluster ID, coordinates, g magnitude, g—r and g—i colours are from our GMOS imaging 
and arc instrumental magnitudes only. Heliocentric velocities arc from the spectra presented in this work. 
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of Thomas, Maraston & Korn (2004; hereafter TMK04) were chosen because they present 
the only set of models that include the effect of a abundance ratios on the Balmer lines. 

We compare the measured Lick indices to the TMK04 SSPs and obtain a minimum x 2 fit. 
This fit is obtained using all the indices measured. Simultaneously a set of \ 2 minimisation 
fits are found with each of the indices omitted. From this set, we select the fit with the 
lowest total x 2 , remove the necessary index and repeat until a stable fit is achieved with no 
highly aberrant (> 3er) indices remaining. All GCs had some indices that were significant 
outliers to the fit and therefore removed during this process. For gl426 none of the indices 
potentially affected by the PN style emission lines are included in the final fit. 

The molecular band indices Mgi and Mg 2 are systematically offset due to poor flux 
calibration and were excluded for all GCs (see Proctor et al. 2005). Similar to other GC 
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Table 2. Globular Cluster indices A < 4600 A. Central index values (first line) and errors (second line). Indices in brackets are removed during 
the fitting process. Missing values are due to limited wavelength coverage. 
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gl078 
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0.588 
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0.195 
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0.305 
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0.356 
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2.917 
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(0.012) 


(0.014) 


(0.221) 
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0.357 
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0.259 


(0.398) 
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(-0.037) 


(-0.008) 


(-0.208) 
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0.051 
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(0.011) 


(0.013) 
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(0.475) 
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-1.286 


6.183 


(-1.965) 


(5.771) 






(0.927) 


(0.025) 


(0.029) 


(0.339) 


0.650 


0.851 


0.487 


0.960 


(0.555) 


(0.615) 
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Table 3. Globular Cluster indices A >4600 A. Central index values (first line) and errors (second line). Indices in brackets are removed during 
the fitting process. Missing values are due to limited wavelength coverage. 



ID 




n P 




IVlgl 


A/Tern 
lvlg2 


Mgb 


Fc5270 


Fp^^ 




Fo^7flQ 


Fc5782 




(A) 


fAl 


(A) 


l md gj 




fAl 


fAl 


(A) 


fAl 


fAl 


fAl 


gl566 


(1.474) 


1.994 


2.692 


(0.039) 


(0.110) 


1.895 


1.668 


(1.606) 


0.832 


0.467 


0.242 




(0.180) 


0.065 


0.139 


(0.001) 


(0.002) 


0.064 


0.070 


(0.080) 


0.059 


0.044 


0.042 


g296 


0.339 


(2.679) 


1.676 


(-0.039) 


(0.028) 


(1.027) 


1.194 


0.819 


(0.778) 








0.382 


(0.132) 


0.283 


(0.003) 


(0.003) 


(0.130) 


0.142 


0.163 


(0.119) 






gl940 


1.193 


2.153 


2.816 


(0.008) 


(0.093) 


2.045 


1.745 


1.425 


(1.250) 


(0.772) 


0.178 




0.362 


0.131 


0.279 


(0.003) 


(0.003) 


0.128 


0.142 


0.161 


(0.117) 


(0.087) 


0.084 


gl550 


(3.430) 


(1.100) 


3.697 


(0.085) 


(0.198) 


2.810 


2.206 


1.687 


1.354 


0.727 


0.391 




(0.369) 


(0.139) 


0.287 


(0.003) 


(0.003) 


0.133 


0.143 


0.163 


0.118 


0.089 


0.085 


g2122 


2.411 


1.670 


3.940 


(0.065) 


(0.200) 


2.953 


2.387 


(2.850) 


1.546 


0.590 


(0.346) 




0.394 


0.144 


0.304 


(0.003) 


(0.004) 


0.139 


0.150 


(0.165) 


0.123 


0.093 


(0.089) 


gl078 


(3.035) 


2.147 


(4.188) 


(0.025) 


(0.112) 


2.118 


1.495 


1.113 


1.115 


0.685 


0.314 




(0.366) 


0.134 


(0.280) 


(0.003) 


(0.003) 


0.130 


0.144 


0.165 


0.119 


0.087 


0.084 


gl420 


1.273 


(0.938) 


(-0.436) 


(0.038) 


(0.115) 


1.418 


1.153 


1.115 


0.694 


0.250 


0.302 




0.380 


(0.143) 


(0.296) 


(0.003) 


(0.003) 


0.136 


0.147 


0.166 


0.121 


0.093 


0.088 


gl540 


(-1.675) 


2.934 


(-0.446) 


(0.016) 


(0.056) 


1.441 


0.717 


0.898 


(0.973) 


(-0.433) 


-0.156 




(0.414) 


0.139 


(0.325) 


(0.003) 


(0.004) 


0.144 


0.163 


0.184 


(0.133) 


(0.101) 


0.095 


gl617 


1.517 


(1.131) 


2.076 


(-0.032) 


(0.076) 


2.534 


1.268 


0.814 


0.354 


0.288 


0.163 




0.448 


(0.169) 


0.352 


(0.003) 


(0.004) 


0.155 


0.176 


0.203 


0.148 


0.112 


0.106 


gl588 


(-0.540) 


2.255 


(-0.892) 


(0.081) 


(0.167) 


2.413 


(3.422) 


1.600 


(-0.020) 


0.700 


(0.002) 




(0.546) 


0.187 


(0.418) 


(0.004) 


(0.005) 


0.185 


(0.195) 


0.226 


(0.175) 


0.121 


(0.118) 


gll4 


(2.781) 


2.327 


2.621 


(0.001) 


(0.051) 


1.230 


1.687 


1.310 


0.706 








(0.558) 


0.197 


0.426 


(0.004) 


(0.005) 


0.196 


0.212 


0.245 


0.180 






gl426 


0.957 


2.962 


(-0.001) 


(-0.067) 


(-0.005) 


2.373 


(-1.293) 


-0.045 


0.920 


(-0.364) 


0.161 




0.570 


0.196 


(0.433) 


(0.004) 


(0.005) 


0.188 


(0.240) 


0.267 


0.188 


(0.143) 


0.131 


gl610 


0.796 


2.451 


3.175 


(0.024) 


(0.075) 


(0.797) 


0.906 


0.900 


0.419 








0.576 


0.207 


0.441 


(0.004) 


(0.005) 


(0.208) 


0.222 


0.252 


0.189 






gl571 


1.771 


0.976 


(0.318) 


(0.085) 


(0.157) 


2.159 


2.677 


2.220 


(0.742) 


(-0.124) 


0.705 




0.578 


0.214 


(0.464) 


(0.004) 


(0.005) 


0.207 


0.221 


0.249 


(0.192) 


(0.138) 


0.126 


gl519 


0.895 


(0.468) 


(6.229) 


(-0.003) 


(0.189) 


(4.168) 


(3.613) 


1.874 


(-0.074) 


(1.240) 


0.362 




0.549 


(0.213) 


(0.406) 


(0.004) 


(0.005) 


(0.180) 


(0.190) 


0.226 


(0.173) 


(0.129) 


0.126 


gl934 


0.258 


3.209 


(-2.122) 


(-0.023) 


(-0.052) 


(-3.180) 


1.031 


0.496 


(0.961) 


0.157 


-0.198 




0.662 


0.231 


(0.533) 


(0.005) 


(0.006) 


(0.260) 


0.264 


0.301 


(0.216) 


0.168 


0.167 


gl580 


2.831 


(3.334) 


(7.071) 


(0.042) 


(0.166) 


(3.601) 


2.713 


(0.901) 


(0.421) 


0.856 


(-0.083) 




0.696 


(0.246) 


(0.518) 


(0.005) 


(0.006) 


(0.236) 


0.260 


(0.306) 


(0.225) 


0.163 


(0.161) 


g219 


(-2.680) 


2.595 


0.108 


(-0.034) 


(0.043) 


1.345 


0.930 


-0.152 


0.715 








(0.809) 


0.274 


0.614 


(0.006) 


(0.007) 


0.271 


0.299 


0.346 


0.246 






gl893 


7.066 


(3.030) 


7.438 


(0.111) 


(0.250) 


5.102 


(4.318) 


2.351 


2.286 


1.414 


(0.316) 




0.818 


(0.305) 


0.615 


(0.006) 


(0.007) 


0.273 


(0.293) 


0.345 


0.254 


0.191 


(0.182) 


g!630 


(1.358) 


2.091 


(0.673) 


(0.145) 


(0.238) 


2.505 


2.254 


2.013 


(-0.329) 


1.568 


1.053 




(0.979) 


0.354 


(0.827) 


(0.007) 


(0.009) 


0.351 


0.366 


0.417 


(0.322) 


0.219 


0.212 
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Table 4. Derived Globular Cluster Properties. Age, [Fe/H], [E/Fe] and [Z/H] arc derived from the x 2 minimisation process, with errors derived 
by a Monte Carlo style method. [Fc/H] ^ H 1S derived according to the method of Brodic & Huchra (1990J from a reduced sample of indices. 



ID 


Age 


[Fe/H] 


[E/Fe] 


[Z/H] 


[Fc/H] BH 




(Gyr) 


(dex) 


(dcx) 


(dex) 


(dcx) 


glDOO 


1 O A4-1 1 
1 Z . zt 1 . 1 


i 1 7-UA 
—1.1 f ztu.uo 


n 1 c:_i_n nc 
U. lOltU.UO 


1 A*34-A HQ 

— l.UoztU.Uo 


1 1 K _|_ A Afl 
— l.lOltU.UD 


gzyo 


1 1 9-1-1 K 


— ± .oZzLU. 1 1 


n 71 4-0 1 n 


1 i f;4_n flfi 

1 . lOltU.UO 


1 SA4-0 9fi 
— 1 .04ICU.ZD 


n 1 CM a 
giy4U 


1 1 Q-LO 1 


1 9fi4-fl 1 1 
— 1 .ZO^U. 1 1 


n ^fi4-n 1 n 
U.Ouztu. lu 


— u.yoznu.uo 


1 Ql 4_A 1 O 
— 1 .OllCU. 1Z 


... 1 Ken 

gioou 


10.U±4.o 


— U.oy±U. 11 


a 1 c^4-a n7 

U. lOltU.U f 


n c;c;4_a 1 1 
— U.OOltU.ll 


A e;04_A 1 fi 
— U.OZltU.lo 


<v91 99 

gzizz 


1 £ A4-9 n 
lo.Uztz.U 


— U.4 / ztu.U / 


a 1 k4_a nfi 

U. lOltU.UO 


A QQ4-A A/1 

-U.Ov5itU.U4 


A A c^4-A 1 A 
— U.40±U.1U 


glU I O 


1 fl 04-1 7 
lU.UZtl. i 


— i .uoztu.uy 


u.z4zcu.uy 


— U.OOICU.UD 


1 1 74-A Afi 
— 1 . 1 I ICU.UD 


... i /I on 
g!4ZU 


1 a o —i— "i a 

lU.Uzt 1. 1 


i /in_i_n 1 1 

— 1 .4UltU. 1 1 


a 91 4-a no 
u.ziztu.uy 


1 OA4-A Afl 
— l.zUltU.UD 


1 OQ J_A 1 A 
— l.ZOltU.lU 


glD4U 


Q A -1-1 7 

y.4± i. ( 


o ao4-a 1 s 
— z.uy±U. lo 


n Afi4-n 1 7 

U.OoltU. 1 ( 


1 /I C^4-A no 

— i.4o±u.uy 


1 fiA-\-C\ C\A 
— l.D4±U.U4 


glol / 


1 1 94-1 o 


— l.ozztu. lz 


n Afi4-n 1 a 

U.OoltU. 1U 


A 9.S4-A Afi 
— U.ooltU.Uo 


1 c;o4-A oc 
— i.oy±u. zo 


n 1 coo 


£ A 4-9 Q 


fl ^74-0 1 A. 
— U.O I ZlU. 14 


n 1 c:4_n HQ 

u. loztu.uy 


fl 4^4-0. 1 ^ 
— U.40ICU. lo 


fl ^34-0 1 
— U.OOlCU. 10 


rr~\~\ A 


lo.Uztu.0 


11/1 _i_n 1 fi 
— 1 . 14ztU. ID 


n AA4-A i c; 
— U.UOltU. ID 


1 OA4-A 1 A 
— 1 .ZUltU. 14 


1 AA4-A 1 7 
— l.OUltU. 1 i 


g!4Z0 


7 c: _|_o c 


i c;Q_i_n 1 q 
— l.OoztU. lo 


fl Q A4- A fl7 
U.oUltU.U / 


A Q Q J_ A 1 1 
— U.oo±:U. 11 


9 Q'74-A AA 
— Z.O / ±U.44 


gl610 


10.6±2.1 


-1.62±0.13 


0.71±0.12 


-0.95±0.09 


-1.98±0.40 


gl571 


15.0±3.6 


-0.10±0.08 


-0.27±0.05 


-0.35±0.08 


-0.29±0.21 


gl519 


12.6±2.9 


-0.72±0.13 


-0.03±0.10 


-0.75±0.09 


-0.82±0.37 


gl934 


10.0±2.3 


-2.31±0.31 


0.65±0.29 


-1.70±0.15 


-2.27±0.21 


gl580 


9.4±2.5 


-0.24±0.13 


0.09±0.09 


-0.15±0.09 


-0.73±0.12 


g219 


12.6±1.8 


-2.72±0.37 


0.50±0.35 


-2.25±0.13 


-2.15±0.25 


gl893 


11.2±5.9 


0.12±0.16 


0.24±0.09 


0.35±0.13 


0.44±0.05 


gl630 


11.2±5.4 


0.10±0.19 


-0.24±0.08 


-0.13±0.19 


0.20±0.35 


gl595 


8.9±3.2 


-3.00±0.49 


0.80±0.41 


-2.25±0.31 


-2.08±0.42 


g230 


15.0±5.9 


-1.52±0.30 


-0.30±0.34 


-1.80±0.23 


-2.34±0.94 



studies (e.g. Beasley et al. 2004a) we find the CN indices to be enhanced relative to the 
models and therefore they were also removed. 

An index- index plot of B.'Ja vs Fe4383 is shown in Fig. 2. This pairing of indices was 
chosen as the combination of an age and a metallicity sensitive index that had the largest 
number of reliable points. This shows the majority of GCs lie near the 15 Gyr age line. 
Horizontal branch morphology is the apparent source of the offset to the left of the TMK04 
15 Gyr age line seen in the majority of the low metallicity GCs (see Pierce et al. 2005 
for a brief discussion on the effect of horizontal branch morphology on GC x 2 fitting). 
TMK04 models in this regime posses a predominantly blue horizontal branch whereas the 
red horizontal branch models from Thomas, Maraston & Bender (2003) pass directly through 
these points. There are a couple of clear outliers for which one or both index values have 
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Figure 2. The indices vs Fe4383 are plotted against agc-mctallicity grids from TMK04. We show the grid for [E/Fe]=+0.3, 
with mctallicity in 0.25 dex steps from from —2.25 (left) to +0.5 (right) and ages of 1,2,3,5,8,12 and 15 Gyrs (top to bottom), 
the bold line is 15 Gyrs and the dashed line is the red horizontal branch 15 Gyr isochrone from Thomas, Maraston & Bender 
(2003). Filled circles are shown for GCs where both indices are included in the fits, if either index is excluded then an open 
circle is plotted. A large number of the GCs fall near the maximum age line for the Thomas, Maraston & Bender (2003) red 
horizontal branch SSPs. 



Figure 3. GC age-metallicity relation. The plot shows that the observed GCs are consistent with an old age and a wide range 
of mctallicities. 



been excluded from the fitting process. However to derive values we use the results from the 
X 2 fitting process. 

From our quantitative \ 2 minimisation fits we find all of the GCs to be consistent with 
an old age ( > 10 Gyrs) with a spread of metallicities from [Fe/H] < -2 to solar (Fig. 3). 
This is consistent with Fig. 2 and indeed the fitting procedure has identified clear outliers. 
The GC, gl426, with the youngest fit age of 7.5 Gyrs, is the GC with unexplained emission 
lines in its spectrum. It also has [E/Fe]=+0.8 which is at the maximum of the models. While 
none of the clearly affected indices were included in the fit, we suspect that the fit was still 
somewhat influenced by the emission source. 

A test of SSP-derived metallicities is to compare them with those derived using the 
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-3 -2 -1 
[Fe/H] 

Figure 4. A comparison of Brodic-Huchra mctallicity vs our derived metallicity. The four points with poor agreement at low 
mctallicity (below [Fc/H]bh= -2) all are at the limits of the SSP models. A gap in the metallicity distribution can be seen at 
[Fc/H]~ -1, which resembles that of the Milky Way GC system. 




Observed g — i 

Figure 5. The predicted colour from stellar population fits vs observed g-i colour. The outliers are GCs that are generally 
at the limits of the SSP models. The errors for both observed and predicted colours arc lower limits and may be significantly 
underestimated. The dashed lines show the RMS scatter of ±0.1 mag around the one-to-one solid line. A Spearman Rank test 
gives a probability of correlation of 99.2%. 

Brodie & Huchra (1990; hereafter BH) method. We were unable to use the full sample 
of metal-sensitive indices due to the poor sensitivity at shorter wavelengths. We therefore 
measured the G band, Fe52, MgH and Mg 2 BH indices using the method outlined in their 
paper. We define a Brodie-Huchra metallicity as the average of the empirically calibrated 
metallicity from these 4 indices, where the error quoted is the standard deviation. This 
metallicity is referred to as [Fe/H] B # in Table 4. Fig. 4 shows that the BH metallicities 
closely match those derived by x 2 fitting to SSPs. A gap in both metallicity measures occurs 
at [Fe/H]~ -1. This is similar to that observed in the Milky Way GC system (Harris 1996), 
which is also composed of an old GC population. The outliers at low metallicity are at the 
limits of the SSP models for either [Fe/H] or [E/Fe]. 
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Figure 6. A plot of a-element ratio vs metallicity. The four GCs with poorly determined mctallicitics and a-element abundance 
ratios (i.e. [E/Fc] errors >0.3) are not shown. The plot indicates a correlation in the sense of increased enhancement with 
decreasing metallicity. 

We compare the observed g-i colours for the GCs with the predicted g-i colours from the 
Bruzual & Chariot (2003) SSPs, with the ages and metallicities from TMK04 SSPs in Fig. 

5 (thus there could be some systematic offsets due to model differences). The use of Bruzual 

6 Chariot (2003) colour models is necessary because the TMK04 SSPs do not include g 
and i magnitudes. There are several outliers in the colour-colour plot, these include the low 
metallicity GCs mentioned previously as well as a few others that are at the SSP model 
maximum age. A Spearman Rank test gives the probability of correlation as 99.2%. Overall 
the figure gives us confidence in the SSP-derived parameters. 

TMK04 estimate a-element abundance ratios using the parameter [E/Fe]. A definition 
of [E/Fe] for the SSP models used can be found in Thomas et al. (2004). Briefly, it includes 
a-elements, such as O, Ne, Mg, Si, S, Ar, Ca and Ti plus two non-a-elements N and Na. 
The [E/Fe] vs. [Fe/H] plot presented in Fig. 6 shows decreasing a-element abundance ratio 
with increasing metallicity, including two solar metallicity GCs with sub-solar a-element 
abundance ratios. The source of this trend is confirmed by the index-index plots of Fig. 7, 
in which high-metallicity GCs generally have a-sensitive indices that are consistent with 
[E/Fe]<+0.3. The lower metallicity GCs have a-element sensitive indices that are predomi- 
nantly stronger than [E/Fe]=+0.3. 

4 GLOBULAR CLUSTER KINEMATICS 

In this section we use our measured velocities for 22 GCs (see Table 1) to estimate the 
radial velocity dispersion and mass-to-light ratio for the halo of NGC 3379. Our sample is 
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Figure 7. The upper row show the a-element sensitive Mgb index vs three different Iron indices (Fc4383, Fe4531 and Fc5270). 
The middle row shows C4668 vs the same Iron indices. All six have TMK04 12 Gyr isochrones plotted with [E/Fe]=+0.3 as 
solid and [E/Fc]=+0.0 as dashed lines. Filled symbols show where both indices are included in the fit, if either of the indices 
are not included in the final fit then open symbols are shown. For several plots one or two highly aberrant indices fall outside 
the range plotted. All six plots show that the «-sensitive indices are above the [E/Fc]=+0.3 dex model at low Iron sensitive 
index values and arc below for high Iron values. The bottom row of plots show the index offset relative to the [E/Fe]=+0.3 dex 
model line, for included indices. Square symbols are Mgb and the circles C4668. A strong trend of decreasing a-element ratio 
with increasing Iron abundance is seen for both o-element sensitive indices relative to the [E/Fe]=+0.3 model expectation. 

supplemented by measurements of 14 GC velocities from Puzia et al. (2004), who used the 
FORS instrument on the VLT. 

The spatial distribution of the combined 36 GCs, in relation to NGC 3379 and NGC 3384, 
is shown in Fig. 8. From this figure, it can be seen that four of the GCs (from our GMOS 
sample) lie roughly half way between NGC 3379 and NGC 3384. As NGC 3379 is more 
luminous than NGC 3384 (My -21.06 vs. -20.64), we might expect NGC 3379 to have 
the more populous GC system and indeed this seems to be the case (Whitlock, Forbes & 
Beasley 2003). The average velocity of the four GCs is 924 ±46 km/s. The systemic velocity 
of NGC 3379 is 911 km/s. For NGC 3384, the systemic velocity is 704 km/s; however, 
NGC 3384 has a substantial rotation (it is an SO galaxy). The rotation amplitude is about 
200 km/s (Fisher 1997) and the side closest to NGC 3379 is redshifted. Any rotation in 
NGC 3379 GC system is below 2 a significance at all radii and the sample is consistent with 
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Figure 8. Spatial positions of GCs. The circles are GMOS data presented here, crosses are VLT data of Puzia et al. (2004). The 
two large circles show the effective radius of the galaxies NGC 3379 (lower right) and NGC 3384 (upper left) (de Vaucouleurs 
et al. 1991, Capaccioli et al. 1990). The effective radius for NGC 3379 is 55" (2.8 kpc). The four GCs approximately midway 
between NGC 3379 and NGC 3384, most likely belong to NGC 3379, however they cannot be decisively assigned to either 
galaxy and are therefore shown as open circles. 



no rotation. In the unlikely situation of the NGC 3384 rotating disk extending to several 
effective radii, and that the four GCs in question are part of the NGC 3384 disk population, 
the expected mean velocity of NGC 3384 members would be essentially the same as that for 
NGC 3379. Therefore, the membership of these four GCs is not entirely clear. We run the 
dynamical models below, both including and excluding these 4 GCs, but our conclusions are 
not greatly affected either way. 

Ideally, the best estimate of the mass-to-light (M / L) variation would come from models 
without assumptions about the distribution function. The state-of-the-art is to use orbit- 
based models (e.g., Gebhardt et al. 2003; Thomas et al. 2005). However, with only 36 clusters 
it is not practical to run these flexible models. For an estimate of the Mj L variation at large 
radii, we instead employ isotropic models. While there are theoretical concerns for assuming 
isotropy at large radii, these models provide a comparative base to other studies. 

For an isotropic model, we only require the second moment of the projected velocity 
distribution. The data on the inner parts of the galaxy come from Gebhardt et al. (2000) 
and Statler & Smecker-Hane (1999). For the velocity dispersion profile from the GC system 
we use a lowess estimator (explained in detail in Gebhardt et al. 1995). The lowess technique 
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R (arcsec) 

Figure 9. Absolute velocity differences between GCs and the systemic velocity of NGC 3379 (911 km/s) with projected 
galactocentric radius. The circles are GMOS and VLT data presented here, open circles indicate the four GCs which are 
possibly associated with NGC 3384. The red triangles are the binned PN data of Romanowsky et al. (2003). The hatched blue 
region shows the 1 a confidence interval for the velocity dispersion from the lowess estimator, including all GCs. 



effectively runs a radial window function through the data to estimate the velocity squared. It 
takes the velocity uncertainties into account directly We have checked this estimator against 
a variety of dispersion estimates, including maximum likelihood, and find no significant 
differences. 

Fig. 9 shows the line-of-sight velocity with projected radius for the individual GCs. This 
allows a comparison of our unsmoothed data with the lowess estimator velocity dispersion 
and the binned Romanowsky et al. (2003) PN data. The individual data point at R~150"with 
| AV|~400 km/s has a relatively large velocity difference which increases the velocity dis- 
persion at large radii. If we remove this one point, then the dispersion is reduced from 166 
km/s to 125 km/s at the largest radii (which are the most uncertain). When compared to 
the PN dispersion, this results in a 3.5-a difference as opposed to a 4-ct difference when the 
full sample is used. The technique we are using for the dispersion has some robustness to 
outliers built in already, which is why the change in statistical significance is not particularly 
strong. 

The top panel in Fig. 10 plots the dispersion profile combining both the stellar and 
globular cluster data. In the overlap region the agreement in the dispersion estimates is 
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very good. The 68% confidence bands come from Monte Carlo simulations as detailed in 
Gebhardt et al. (1995). 

We use non-parametric, isotropic models as outlined in Gebhardt & Fisher (1995). Given 
the surface brightness profile and velocity dispersion profile, the spherical Jeans equation can 
uniquely determine the mass density profile (and hence the Mj L profile) assuming isotropy. 
Through the Abel deprojection, the surface brightness profile uniquely determines the lu- 
minosity density. Similarly, the surface brightness times the projected velocity dispersion 
determines the luminosity density times the internal (3-D) velocity dispersion. From the 
internal velocity dispersion, we can derive the mass profile and, hence the M/L profile. We 
employ a degree of smoothing since the Abel deprojection involves a derivative (Gebhardt 
et al. 1996). 

The surface brightness profile for the kinematic tracer is important to characterize well 
for the dynamical analysis. We use the globular cluster number density profile from Rhode & 
Zepf (2004). Their data extend from 1.2' to approximately 20'. The surface brightness profile 
for the stellar light comes from Gebhardt et al. (2000) which extends to 2.8'. Thus, there is 
significant overlap for comparison. In the overlap region, the stellar light is slightly steeper 
than the GC profile, and after 2.5' the GC profile flattens significantly. For the dynamical 
analysis, one should use the GC profile only; however, in order to de-project properly, one 
requires the central profile which does not exist for the GCs. Thus, we have to rely somewhat 
on the stellar profile for the extrapolation inward. For the dynamical analysis, we use the 
GC profile from Rhode & Zepf into 1.2', and then use the stellar profile interior to that. 
We have also tried a variety of profiles: using only the cluster profile with various central 
extrapolations, using only the stellar light profile, and using various extrapolations to larger 
radii. We find that there is essentially no effect from the large radii extrapolation. For the 
small radii extrapolation, the effect on the projected dispersion is dramatic at small radius 
but relatively insignificant at the radii where we have GC velocities. Thus, all of our tests 
suggest there is little effect on our overall results from the surface brightness profile. 

The velocity dispersion and M/L ratio are plotted with galactocentric radius in Fig. 10. 
GC data, both including and excluding the four GCs potentially associated with NGC 3384, 
are shown. The heavy solid line in the upper panel in Fig. 10 is the expected isotropic velocity 
dispersion profile for a constant M/L ratio, which is consistent with the Romanowsky et al. 
(2003) PN data. In the bottom panel we plot the projected M/L profile, calculated by the 
lowess method and using the surface brightness profile characterization described above. 
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Figure 10. Plots of radial velocity dispersion and M/L ratio using GC and stellar data (from Gebhardt et al. 2000 and Statlcr 
& Smecker-Hane 1999). The upper panel shows the velocity dispersion of the GCs. Blue open squares and hashing include the 
four GCs that that may be associated with NGC 3384, and the red hashing and filled squares when they are excluded. The 
squares show the binned data and the hashed areas show the range of values using the lowess estimator (with 1 a errors). The 
PN data points from Romanowsky et al. (2003) (circles) and the expected isotropic profile with constant M/L (thick line) are 
also shown. In the lower panel, the blue line (1 a hatching) shows that when the four ambiguous GCs are included the M/L 
ratio rises slowly. The red line (with 1 a hatching) shows that the M/L ratio rises strongly in the outer regions when the four 
ambiguous GCs are excluded. The horizontal line in the bottom panel is M/L = 6. See text for details. A colour version of 
this figure is available in the online journal. 

Both cases (including and excluding the four GCs) shows a rise in the M/L at large radii. 
This rise in M/L is stronger when the four GCs are excluded. 

Finally, we remind the reader that the above results are based on only 36 GCs as kine- 
matic tracers; clearly such an analysis should be repeated when larger samples become 
available. 

5 DISCUSSION 

The reliability of GC parameters derived from Lick indices and our \ 2 fitting method are 
affected by several factors, most importantly the stability of the fit when individual indices 
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are excluded. For some GCs it is not possible to find a stable fit that includes a large number 
of indices (>10). This is often the case for the low S/N GCs. The large index errors of low 
S/N GCs mean that, within 3 a, the indices can be consistent with a wide range of ages and 
metallicities. In this situation the fit can be driven by a single index which is outside the 
range spanned by the models, leading to a fit at the extreme of the SSP models. When this 
occurs the Monte Carlo determined errors are large, reflecting the lower accuracy of these 
parameters. 

To examine the effects of these factors we tested the reliability of the x 2 At metallicities 
by comparing them to Brodie-Huchra metallicities. Fig. 4 shows the good agreement between 
the two metallicity measures, with some divergence below [Fe/H]=-2 dex. We then used g-i 
colour as a test of the derived age and metallicity. Model predictions for g-i from BC03, 
for our x 2 At ages and metallicities (derived using TMK04 SSPs), are plotted against the 
observed g-i colours in Fig. 5. We see reasonable agreement between the predicted and 
observed colours, with the outliers being those low S/N GCs that are at the limits of the 
models. 

Having established the reliability of our x 2 minimisation fits for age and metallicity, we 
find, all of the GCs are consistent with an old age ( > 10 Gyrs) within errors and a spread 
of metallicities is seen from [Fe/H] <-2 to solar (see Fig. 3). We see no discernible age 
structure. 

We next consider the a-element abundance ratio parameter [E/Fe]. Due to the necessary 
exclusion of the Mg and CN indices (see Section 3), the only strongly a-element sensitive 
indices included in the fits are Mgb and C4668. However, all the plots of Mgb and C4668 
against Fe sensitive indices (Fig. 7) display the same overall trend, of decreasing a-element 
ratio with increasing metallicity. This is highlighted by plotting the residuals of the model 
[E/Fe]=+0.3 line to these index pairings. This is also seen for [E/Fe] derived from the 
X 2 fitting process (Fig. 6). A similar trend is apparent for GCs in NGC 1052 (See Fig. 
9 of Pierce et al. 2005) and is also noted in Puzia (et al. 2005). This contrasts with the 
findings of Carney (1996) for Galactic GCs in which [ce/Fe] is largely constant for GCs of all 
metallicities. However, the Puzia et al. (2005), NGC 3379 and NGC 1052 results are derived 
from TMK04 models. These models were constructed using the results of Trippico & Bell 
(1995) which included the effects of a-element variation on stars of solar-metallicity. The 
application of the Trippico & Bell (1995) models to GCs of [Fe/H]~-2 may therefore play 
a part in generating the observed trend of high a-element ratios for GCs at low metallicity. 
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At such low metallicity we can see from Fig. 7 that the indices cannot offer much leverage 
to differentiate the [E/Fe] values. 

However, we note that the a-element ratio trend extends up to the solar metallicity at 
which the Trippico & Bell (1995) models were calculated. Of particular interest are the two 
GCs with solar metallicity, [E/Fe]~-0.25 and ages greater than 10 Gyrs. While there is no 
obvious link, it is interesting to note that X-ray observations of hot gas in the centre of the 
NGC 5044 galaxy group reveal a similar abundance pattern, of a sub-solar a-element ratio 
at solar metallicity (Buote et al. 2003). A significant number of type la, rather than type 
II, supernovae are therefore necessary to produce enough Iron to explain this abundance 
pattern. Such a process must have occurred rapidly due to the old ages of the GCs. 

From the measured recession velocities of our sample of 22 GCs, and 14 GCs from Puzia 
et al. (2004), we find evidence for a significant dark matter halo. This result does not depend 
on the inclusion of 4 GCs which may belong to NGC 3384. The PN data of Romanowsky 
et al. (2003) showed a decrease in velocity dispersion at large radii, suggesting a minimal 
DM halo. However, Dekel et al. (2005) recently showed that stellar orbits in the outer 
regions of merger-remnant elliptical galaxies are elongated and that declining PN velocity 
dispersions do not necessarily imply a dearth of dark matter. This demonstrates that the 
orbital properties of the kinematic tracers need to be well understood. For PN radial orbits 
need to be included, while for GCs isotropic models should be sufficient (see Cote et al. 
2003), but caution is still required. Ignoring the differences between dynamical models, we 
still find the results from the two studies contradict each other. 

A two-sided x 2 test over an radii using all GCs, gives a 2 a difference between the mass 
profile derived from our GCs compared to the PN-derived mass model of Romanowsky et 
al. (2003). This would indicate borderline statistical difference between our mass model and 
the no-dark matter model of Romanowsky et al. . 

To test the significance of the differences between the PN and GC data we have carried 
out a Monte Carlo simulation. At the radius of a given PN binned data point, we sample 
from a Gaussian distribution using a mean given by the GC Lowess dispersion estimate 
and the standard deviation from the 68% confidence band at that point. We draw another 
Gaussian from the PN value and its uncertainty. We then ask whether the new PN dispersion 
is smaller than the new GC dispersion at that point. For a single measurement, this is only 
a x 2 test; however, generalizing for a range of radii, we determine whether each of the PN 
dispersions are smaller than each of the GC dispersions. This test is rather strict in the 
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sense that if any realisation shows one PN value above the GC value at that radius then 
the difference is not considered. At small radii, the PN and GC dispersions are consistent, 
whereas at large radii, the PN dispersions are lower. 

We find that the PN dispersions are smaller than the GC dispersions at >97% confidence 
beyond an effective radius. Since the PN dispersions are consistent with mass following light, 
the GC dispersions suggest a need for a dark halo. Either the orbital properties of the two 
tracers are different (Dekel et al. 2005), or one or both are hampered by low number statistics. 

Two alternate observations suggest the presence of a dark matter halo around NGC 3379. 
Chandra archival data for NGC 3379 indicate that a hot gaseous halo is present (Fukazawa et 
al. 2006). The kinematics of the HI gas ring around the Leo Triplet of NGC 3379, NGC 3384 
and NGC 3389 suggests a mass-to-light ratio of 27 which is consistent with a DM halo 
(Schneider 1989). 

One possible explanation for the difference between the PN data and the indication of a 
normal DM halo, could be that NGC 3379 is a face-on SO galaxy, as suggested by Capaccioli 
et al. (1991). If a significant fraction of the PN belong to a disk, this would suppress the 
line-of-sight velocity dispersion of the PN relative to that of the GCs which lie in a more 
spherical halo. 

6 CONCLUSIONS 

We have obtained Gemini/GMOS spectra for 22 GCs around the elliptical galaxy NGC 3379. 
We present ages, metallicities and a-abundance ratios that were derived by applying the 
multi-index \ 2 minimisation method of Proctor & Sansom (2002) to the SSP models of 
Thomas et al. (2004). Metallicity estimates, derived according to the method of Brodie & 
Huchra (1990), agree closely with those from our \ 2 minimisation method. We also find 
good agreement between the observed colours and those predicted from our x 2 minimisation 
method ages and metallicities. All the GCs are found to be consistent with old ages, i.e. 
> 10 Gyr, with a wide range of metallicities. This is consistent with the resolved stellar 
population work of Gregg et al. (2004) who found the galaxy stars in the outer regions to 
be old. We find no evidence for a young GC sub-population. 

The a-abundance ratios appear to decrease with increasing metallicity, however inter- 
pretation of this trend is complex and requires further work. 

Using the recession velocities of our 22 GCs and 14 GCs from Puzia et al. (2004), we 
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measure the projected velocity dispersion of the GC system and find that it is consistent 
with being constant with radius in the outer regions. With this velocity dispersion profile, 
NGC 3379 appears to posses a dark halo, although we cannot rigorously determine the dark 
halo mass. This is in contrast to the earlier claims of Romanowsky et al. (2003) of a minimal 
DM halo, based on planetary nebulae kinematics. 

A two-sided % 2 test over all radii using all GCs, gives a 2a difference between the mass 
profile derived from our GCs compared to the PN-derived mass model of Romanowsky et 
al. (2003). This would indicate borderline statistical difference between our mass model and 
the no-dark matter model of Romanowsky et al. However, if we restrict our analysis to radii 
beyond one effective radius and test if the GC velocity dispersion is consistently higher, then 
we determine a > 3a difference between the mass models, and hence favor the conclusion 
that NGC 3379 does indeed have dark matter at large radii in its halo. 

7 NOTES ADDED IN PROOF 

We emphasise that the mass model we compare to, and find discrepancy with, is an isotropic 
constant M/L model, which is the simplest interpretation of the PN data implying no DM. 
This is not the same as the preferred mass model presented in Romanowsky et al. (2003), 
which was an orbit model including some dark matter. For detailed comparison to such a 
model, the dynamical characteristics and projection effects for both the GCs and PN would 
need to be taken into account, which is beyond the scope of this work. 

Our attention was drawn to the paper by Samurovic & Danziger (2005). Their analysis 
of the X-ray halo of NGC 3379 leads to predictions of both M/L and velocity dispersion 
with radius that are consistent with those from our GC analysis For example, see their Figs. 
10 and 13 and their statement "We note that beyond 120" (2.2 R e ) a discrepancy between 
PNe estimates and X-ray estimate occurs" . 
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